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ITERATIVE CHANNEL AND INTERFERENCE ESTIMATION AND 

DECODING 

BACKGROUND 

I. Field 

[1001] The present invention relates generally to communication, and more 
specifically to techniques for recovering data in the presence of noise and interference at 
a receiver in a wireless communication system. 

II. Background 

[1002] In a wireless communication system, a transmitter typically encodes, 
interleaves, and modulates (i.e., symbol maps) traffic data to obtain data symbols, which 
are modulation symbols for data. For a coherent system, the transmitter multiplexes 
pilot symbols with the data symbols, processes the multiplexed pilot and data symbols 
to generate a modulated signal, and transmits the signal via a wireless channel. The 
channel distorts the transmitted signal with a channel response and further degrades the 
signal with noise and interference. 

[1003] A receiver receives the transmitted signal and processes the received signal 
to obtain received symbols. For a coherent system, the receiver typically estimates the 
channel response with the received pilot symbols and performs coherent demodulation/ 
detection of the received data symbols with the channel response estimates to obtain 
recovered data symbols, which are estimates of the data symbols transmitted by the 
transmitter. The receiver then symbol demaps, deinterleaves, and decodes the recovered 
data symbols to obtain decoded data, which is an estimate of the traffic data sent by the 
transmitter. 

[1004] In a typical coherent wireless system, the receiver processes the received 
pilot symbols once to obtain the channel response estimates and also performs coherent 
demodulation once on the received data symbols to obtain the recovered data symbols. 
The receiver then performs symbol demapping, deinterleaving, and decoding on the 
recovered symbols in accordance with the coding and modulation schemes used for the 
traffic data. The noise and interference degrade the quality of the recovered data 
symbols and affect the reliability of the decoded data. There is therefore a need in the 
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art for techniques to recover data in the presence of noise and interference at the 
receiver in a wireless communication system. 

SUMMARY 

[1005] Iterative receiver processing techniques that can account for interference and 
provide improved performance are provided herein. These techniques may be used for 
various wireless communication systems and may be implemented in various manners. 
[1006] For an iterative channel and interference estimation and decoding scheme, 
channel gain and interference are iteratively estimated. For this scheme, prior 
information for channel gain and interference is initially obtained (e.g., based on 
received pilot symbols). Forward information for code bits corresponding to received 
data symbols is derived based on the received data symbols and the prior information. 
(Each data symbol is obtained based on B code bits, where B is dependent on the 
modulation scheme used for the data symbol.) The forward information is then decoded 
to obtain feedback information for the code bits corresponding to the received data 
symbols. A posteriori information for channel gain and interference for each received 
data symbol is derived based on the feedback information for the code bits 
corresponding to that received data symbol. The a posteriori information for all 
received data symbols and the prior information are combined to obtain updated 
information for channel gain and interference for each received data symbol. The 
process can be repeated for any number of iterations, with the updated information for 
channel gain and interference being used to derive the forward information for each 
subsequent iteration. 

[1007] The prior information, a posteriori information, and updated information for 
channel gain and interference may be represented with joint probability distributions on 
channel gain and interference. The forward and feedback information may be 
represented by log-likelihood ratios (LLRs) for the code bits of the received data 
symbols. 

[1008] In another iterative receiver processing scheme, interference is estimated 
once and used in the iterative channel estimation and decoding process. This can reduce 
computation complexity for the various steps of the iterative process. 
[1009] For all iterative receiver processing schemes, complexity can be reduced by 
estimating channel magnitude non-iteratively (e.g., based on the received pilot and data 
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symbols) and only estimating channel phase iteratively. To further reduce complexity 
when using an M-ary phase shift keying (M-PSK) modulation scheme, the channel 
phase may be estimated to within a range of 0 to 2niy[ (e.g., based on the received 
data symbols). In this case, the prior information, a posteriori information, and updated 
information for channel gain can comprise M components for M different channel phase 
values, where M may be much fewer than the number of components needed without 
the initial channel phase estimate. 

[1010] Various aspects and embodiments of the invention are described in further 
detail below. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[1011] The features and nature of the present invention will become more apparent 
from the detailed description set forth below when taken in conjunction with the 
drawings in which like reference characters identify correspondingly throughout and 
wherein: 

[1012] FIG. 1 shows a transmission scheme for a frequency hopping system; 

[1013] FIG. 2 shows a transmitter in the frequency hopping system; 

[1014] FIG. 3 shows a receiver in the frequency hopping system; 

[1015] FIG. 4 shows a Tanner graph that graphically illustrates iterative channel and 

interference estimation and decoding; 

[1016] FIG. 5 shows a process for performing iterative channel and interference 
estimation and decoding; 

[1017] FIG. 6 shows a process for performing iterative channel estimation and 
decoding with an interference estimate; 

[1018] FIGS. 7 A and 7B illustrate two methods for estimating channel phase based 
on received data symbols; and 

[1019] FIG. 8 shows a process for performing iterative channel estimation and 
decoding with a channel phase estimate obtained from received data symbols. 

DETAILED DESCRIPTION 
[1020] The word "exemplary" is used herein to mean "serving as an example, 
instance, or illustration." Any embodiment or design described herein as "exemplary" 
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is not necessarily to be construed as preferred or advantageous over other embodiments 
or designs. 

[1021] The iterative receiver processing techniques described herein may be used 
for various wireless communication systems that experience interference. For clarity, 
these techniques are described for a frequency hopping conmiunication system in which 
data is transmitted on different frequency subbands in different time intervals, which are 
also referred to as "hop periods". With frequency hopping, a data transmission hops 
from subband to subband in a pseudo-random manner. This hopping provides 
frequency diversity and allows the data transmission to better withstand deleterious path 
effects such as narrow-band interference, janmiing, fading, and so on. 
[1022] The subbands in a frequency hopping system may be provided by orthogonal 
frequency division multiplexing (OFDM), other multi-carrier modulation techniques, or 
some other constructs. OFDM is a modulation technique that effectively partitions the 
overall system bandwidth into multiple (Np) orthogonal subbands. Each subband is 
associated with a respective subcarrier that may be modulated with data. The subbands 
are also commonly referred to as tones, subcarriers, bins, and frequency channels 
[1023] A frequency hopping system may be deployed with multiple cells, where a 
cell typically refers to a base station and/or its coverage area. Each cell can support 
multiple users simultaneously. For a given cell, the data for each user in the cell may be 
transmitted using a specific frequency hopping (FH) sequence assigned to the user. The 
FH sequence indicates the specific subband to use for data transmission in each hop 
period. Multiple data transmissions for multiple users may be sent simultaneously using 
different FH sequences. These FH sequences are defined to be orthogonal to one 
another so that only one data transmission uses each subband in each hop period. By 
using orthogonal FH sequences, the data transmissions for multiple users in the same 
cell do not interfere with one another while enjoying the benefits of frequency diversity. 
However, these users typically experience inter-cell interference from users in other 
cells. The interference observed by a given user can vary from hop to hop because 
different interfering users may be observed in different hops. 

[1024] FIG. 1 shows an exemplary transmission scheme 100 for a frequency 
hopping communication system. FIG. 1 shows pilot and data transmission on a 
frequency-time plane in which the vertical axis represents frequency and the horizontal 
axis represents time. For this example, Np = 8 and the eight subbands are assigned 



5 



Attorney Docket No. 030372 



indices of 1 through 8. Up to eight traffic channels may be defined whereby each traffic 
channel uses one of the eight subbands in each hop period. A hop period is the time 
duration spent on a given subband and may be defined to be equal to the duration of Nh 
OFDM symbols, where N„ > 1 , 

[1025] Each traffic channel is associated with a different FH sequence. The FH 
sequences for all traffic channels may be generated with an FH function fik^T) , where 
k denotes the traffic channel number and T denotes system time, which is given in units 
of hop periods. Np different FH sequences may be generated with Np different values of 
k for the FH function f(k,T) . The FH sequence for each traffic channel indicates the 
particular subband to use for that traffic channel in each hop period. For clarity, FIG. 1 
shows the subbands used for one traffic channel. This traffic channel hops from 
subband to subband in a pseudo-random manner determined by its FH sequence. 
[1026] For transmission scheme 100, Np pilot symbols (depicted as a solid box) are 
transmitted in a time division multiplexed (TDM) manner with No data symbols 
(depicted as hashed boxes) in each hop period. In general, Np>l, Nj5>l, and 
Np + = Nj^ . Np is typically a sufficient number of pilot symbols to allow a receiver 
to adequately estimate the channel response in each hop period. 

[1027] FIG. 2 shows a block diagram of a transmitter 200 in the frequency hopping 
system. A transmit (TX) data processor 220 receives traffic data from a data source 210 
and control data from a controller 250. Within TX data processor 220, an encoder 222 
encodes the traffic and control data in accordance with a selected coding scheme (e.g., a 
convolutional code, a low density parity check (LDPC) code, a Turbo code, a block 
code, and so on) to obtain coded data. The encoding increases the reliability of the data 
transmission. A channel interleaver 224 interleaves (i.e., reorders) the coded data to 
obtain interleaved data. The interleaving provides diversity for the coded data. A 
symbol mapping unit 226 then symbol maps (i.e., modulates) the interleaved data in 
accordance with a selected modulation scheme to obtain data symbols. The selected 
modulation scheme may be M-PSK (e.g., BPSK or QPSK), M-ary quadrature amplitude 
modulation (M-QAM), or some other modulation scheme. The symbol mapping may 
be performed by (1) grouping sets of B interleaved bits to form B-bit binary values, 
where B > 1 and 2^ = M , and (2) mapping each B-bit binary value to a point in a 
signal constellation corresponding to the selected modulation scheme. Each mapped 
signal point is a complex value and corresponds to a modulation symbol (i.e., a data 
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symbol). Symbol mapping unit 226 provides a stream of data symbols to an OFDM 
modulator 230. 

[1028] OFDM modulator 230 performs frequency hopping and OPDM modulation 
for the data and pilot symbols. Within OFDM modulator 230, a TX FH processor 232 
receives the data and pilot symbols and provides these symbols on the proper subband 
(and in the proper order) in each hop period, as indicated by an FH sequence from 
controller 250. The data and pilot symbols dynamically hop from subband to subband 
in a pseudo-random manner determined by the FH sequence. TX FH processor 232 
provides Nf "transmit" symbols for the Np subbands for each OFDM symbol period. 
These Np transmit symbols are composed of one data/pilot symbol for each subband 
used for data/pilot transmission and a signal value of zero for each subband not used for 
data/pilot transmission. 

[1029] An inverse fast Fourier transform (IFFT) unit 234 receives the Np transmit 
symbols for each OFDM symbol period, performs an Np-point inverse fast Fourier 
transform on the Np transmit symbols, and provides a corresponding "transformed" 
symbol that contains Np time-domain chips. Each chip is a complex value to be 
transmitted in one chip period, where the chip rate is typically determined by the system 
bandwidth. A cyclic prefix generator 236 receives the Np chips for each transformed 
symbol and repeats a portion of the transformed symbol to form an OFDM symbol that 
contains Np + N^p chips, where Ncp is the number of chips being repeated. The 

repeated portion is often referred to as a cyclic prefix and is used to combat inter- 
symbol interference (IS I) caused by a dispersive wireless channel (i.e., a wireless 
channel with time delay spread). An OFDM symbol period is the duration of one 
OFDM symbol, which is Np + N^p chip periods. Cyclic prefix generator 236 provides a 

stream of OFDM symbols. A transmitter unit (TMTR) 242 conditions (e.g., converts to 
analog signals, filters, amplifies, and frequency upconverts) the stream of OFDM 
symbols to generate a modulated signal, which is transmitted from an antenna 244. 
[1030] FIG. 3 shows a block diagram of a receiver 300 in the frequency hopping 
system. An antenna 312 receives the modulated signal transmitted by transmitter 200 
and provides the received signal to a receiver unit (RCVR) 314. Receiver unit 314 
conditions (e.g., frequency downconverts, filters, and amplifies) the received signal and 
further digitizes the conditioned signal to obtain a stream of samples, which is provided 
to an OFDM demodulator 320. 



7 



Attorney Docket No. 030372 



[1031] Within OFDM demodulator 320, a cyclic prefix removal unit 322 receives 
the stream of samples, removes the cyclic prefix appended to each received OFDM 
symbol, and provides a corresponding received transformed symbol that contains Nf 
samples. An FFT unit 324 performs an Np-point FFT on the Nf samples for each 
received transformed symbol to obtain Nf received symbols for the Np subbands for that 
transformed symbol. An RX FH processor/demultiplexer 326 obtains the Np received 
symbols for each OFDM symbol period and provides the received symbol from the 
proper subband as the received data/pilot symbol for that OFDM symbol period. The 
proper subband is determined by an FH sequence from a controller 350. The FH 
sequence used for RX FH processor 326 at receiver 300 is the same as, and 
synchronized to, the FH sequence used by TX FH processor 232 at transmitter 200. RX 
FH processor 326 operates in unison with TX FH processor 232 and provides a stream 
of received data/pilot symbols from the proper subbands to a receive (RX) data 
processor 330. 

[1032] RX data processor 330 performs iterative receiver processing on the received 
data and pilot symbols to obtain decoded data. For the embodiment shown in FIG. 3, 
RX data processor 330 includes a channel and interference estimator 332, a detector 
334, a channel deinterleaver 336, a decoder 340, and a channel interleaver 342 that 
operate as described below. RX data processor 330 provides the decoded data to a data 
sink 348 and/or controller 350. 

[1033] Controllers 250 and 350 direct operation at transmitter 200 and receiver 300, 
respectively. Memory units 252 and 352 provide storage for program codes and data 
used by controllers 250 and 350, respectively. 

[1034] The model for the frequency hopping system may be expressed as: 
r^(m) = Sf^(m)\(m)-^n^{m) 

Eq(l) 

where 5^ (m) is the data or pilot symbol transmitted on subband k in symbol period m; 

\ (m) is the complex channel gain for subband k in symbol period m, which can 
be decomposed into a channel magnitude Uf^im) and a channel phase 

ii(/n) is the received data or pilot symbol on subband k in symbol period m; and 
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rif^ (m) is the noise and interference received on subband k in symbol period m. 

[1035] For simplicity, the channel magnitude (m) = | hf^ (m) | is assumed to be 

known by the receiver and only the channel phase Of^im) needs to be estimated. For a 

constant energy modulation scheme such as M-PSK, the magnitude of data and pilot 
symbols received in each hop period can be averaged to obtain a reasonably accurate 
estimate of the channel magnitude {m) for that hop period. Thus, in the following 

description, the channel gain can be sufficiently characterized by just the channel phase. 
(However, the diagrams show the more general case in which the channel gain can be a 
complex value with unknown magnitude and phase.) 

[1036] The receiver may utilize non-iterative or iterative receiver processing to 
recover the transmitted data. For a non-iterative scheme, the channel response is 
estimated based on the received pilot symbols, and the received data symbols (m) are 
coherently demodulated or "detected" with the channel response estimate to obtain 
recovered data symbols Si^(m)^ which are estimates of the transmitted data symbols 

Sf^ (m) . The detection is performed once for the non-iterative scheme. The recovered 

data symbols are then deinterleaved and decoded to obtain decoded data. For an 
iterative scheme, the channel estimation, detection, and decoding are performed for 
multiple iterations. The iterative scheme exploits the error correction capabilities of the 
coding scheme to provide improved performance. This is achieved by iteratively 
passing information between a channel estimator, a detector, and a decoder for multiple 
iterations, as described below. 

[1037] FIG. 4 shows a Tanner graph 400 that graphically illustrates an iterative 
channel and interference estimation and decoding scheme. Iterative receiver processing 
is performed on a block of data symbols that, in general, can contain any number of data 
symbols. For clarity, the iterative receiver processing is described below for a block of 
Nd received data symbols for one hop period. The No received data symbols are 
formed by Nb code bits, where = B • . 

[1038] Tanner graph 400 includes a channel and interference estimate node 410, Nd 
detection nodes 420a through 420n for Nd data symbols in the block, channel 
deinterleaver 336, channel interleaver 342, decoder 340, and Nd estimation nodes 440a 
through 440n for the Nd data symbols. Node 410 couples to each detection node 420 
via a respective link 412 and to each estimation node 440 via a respective link 442. 
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Each of links 412 and 442 carries information regarding the channel gain and 
interference for the data symbol associated with that link. Each detection node 420 
couples to channel deinterleaver 336 via B links 422, and each estimation node 440 
couples to channel interleaver 342 via B links 436. Decoder 340 couples to channel 
deinterleaver 336 via Nb links 426 and to channel interleaver 342 via Nb links 432. 
Each of links 422, 426, 432 and 436 carries information (in the direction indicated by 
the link) for a code bit associated with that link. The information for each code bit is 
typically in the form of a log-likelihood ratio (LLR), which indicates the likelihood of 
the code bit being a one ("1") or a zero ("0"). 

[1039] For the first iteration, node 410 obtains estimates of the channel gain based 
on the received pilot symbols. Node 410 estimates only the channel phase if the 
channel magnitude is known. Node 410 also estimates the interference observed by the 
received data symbols based on the received pilot symbols. Node 410 provides "prior" 
information for channel phase and interference to each detection node 420 via link 412. 
[1040] Each detection node 420 obtains a respective received data symbol r^(w) and 

the prior information for channel phase and interference from node 410. Each detection 
node 420 computes the IXR for each of the B code bits that forms its data symbol based 
on the received data symbol r^(n) and the prior information for channel phase and 

interference. Each detection node 420 provides B "forward" LLRs for the B code bits 
to channel deinterleaver 336 via links 422. Channel deinterleaver 336 deinterleaves the 
forward LLRs for all Nb code bits in a manner complementary to the interleaving 
performed at the transmitter and provides deinterleaved forward LLRs to decoder 340. 
[1041] Decoder 340 decodes the deinterleaved forward LLRs for the Nb code bits in 
accordance with the coding scheme used by the transmitter. For example, decoder 340 
may implement (1) a maximum a posteriori (MAP) algorithm or a soft-output Viterbi 
(SOV) algorithm if the transmitter uses a convolutional code or (2) a Turbo/LDPC 
decoder if the transmitter uses a Turbo or an UDPC code. Decoder 340 provides 
feedback LLRs for the Nb code bits, which are updated LLRs for these bits, to channel 
interleaver 342 via links 432. Channel interleaver 342 interleaves the feedback LLRs in 
) the same manner as the interleaving performed at the transmitter and provides 
interleaved feedback LLRs to estimation nodes 440 via links 436. 

[1042] Each estimation node 440 obtains a respective received data symbol r^(n) 
and the interleaved feedback LLRs for the B code bits of that received data symbol from 
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channel interleaver 342. Each estimation node 440 derives "a posteriori' information 
for channel phase and interference for its data symbol based on the received data 
symbol (n) and the feedback LLRs and provides this a posteriori information to node 
410 via link 442. 

[1043] Node 410 combines the prior information for channel gain and interference 
and the a posteriori information for channel phase and interference from estimation 
nodes 440 to obtain updated information for channel gain and interference for each 
received data symbol. Node 410 provides the updated information for channel phase 

and interference to each detection node 420, The detection and decoding for the second 
iteration are then performed in similar manner as for the first iteration, albeit with the 
updated information for channel gain and interference. 

[1044] In FIG. 4, each detection node 420 derives and provides "forward" 
information for the code bits corresponding to a respective received data symbol, and 
decoder 340 derives and provides "feedback" information for each received data 
symbol. The forward and feedback information is typically given in the form of LLRs, 
but may also be given in other forms. 

[1045] FIG. 5 shows a flow diagram of a process 500 for performing iterative 
channel and interference estimation and decoding. 

[1046] In step 510, prior information for channel phase and interference is obtained 
based on the received pilot symbols. Step 510 is performed by node 410 in FIG. 4. 
Again, the channel magnitude is assumed to be known and only the channel phase is 
estimated. The normalized channel observation for the pilot symbols may be expressed 
as: 

h = a e'^ + n , Eq (2) 

where the channel gain a and channel phase ^are assumed to be constant for all Nd data 
symbols in the block and the noise and interference n is assumed to be a complex 
Gaussian random variable with zero mean and a variance of Nq . 

[1047] The prior information for channel phase and interference can be given as a 
joint probability distribution on channel phase 0 and interference power /. For 
simplicity, the channel phase 0 may be quantized to L possible values and given in 
integer multiple of Itt/L. Similarly, the interference power / may be quantized to Q 
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possible values. The joint probability distribution on ^and / may be viewed as a three- 
dimensional (3-D) graph with the x-axis representing the channel phase 6^ the y-axis 
representing the interference power /, and the z-axis representing the joint probability of 
a particular channel phase 0^ and a particular interference power given a particular 

received pilot symbol. The joint probability distribution on 0 and / may be expressed 
as: 



where p is the received pilot symbol and p is the actual pilot symbol; 

Pr = ^0,/ = /q I p) is the joint probability distribution on 0 and /, which gives 

the probability of the channel phase being 6^ and the interference power 

being given the received pilot symbol p ; 
Pr = ^0,/ = /q) is an a priori joint probability distribution on 6 and /, which 

gives the probability of the channel phase being 6^ and the interference 

power being 7^; 

Pr(p 1^ = ^0,7 = /(J ) is a probability distribution (obtained based on a 
communication channel model, e.g., Gaussian) that gives the probability 
of obtaining the received pilot symbol p for a given channel phase 6^ 

and interference power 7q ; 
Pr (p) is the probability of obtaining a given value of p ; and 
Pr(7 = 7Q) is an a priori probability distribution on 7, which gives the 

probability of the interference power being . 

[1048] The first expression in equation set (3) is obtained based on Bayes' rule. 
The second expression in equation set (3) is obtained with the assumption that different 
values of channel phase ^are a priori equiprobable, so that Pr (^ = ^o) ^ constant 



Pr(^ = ^o»^ = ^olP) = 



Pr(ff = go,7 = 7o)Pr(p 1^ = ^0,7 = 7o) 
Pr(p) 



oc 



Pr(/ = /o)Pr(p|^ = ^„./ = /o) 



Eq(3) 
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and can be omitted. The third expression in equation (3) is obtained based on the 
assumption that the noise and interference is a complex Gaussian random variable with 
zero mean and variance of 1^ . Equation set (3) omits a normalization factor that can be 

computed by applying the constraint that a probability distribution function (pdf) 
integrates to one over its domain. 

[1049] One joint probability distribution on 0 and / is obtained for each received 
pilot symbol, as shown in equation set (3). If multiple received pilot symbols are 
available, then multiple joint probability distributions on 0 and / are obtained for these 
symbols and combined to obtain one overall or composite joint probability distribution 
on 0 and / for all received pilot symbols. Step 510 provides one joint probability 
distribution on 0 and / for use by all No received data symbols. This joint probability 
distribution contains L Q probability values for L different channel phase values and Q 
different interference power values, which may be viewed as a 3-D graph for probability 
versus 0 and /. This joint probability distribution represents prior information for 
channel phase and interference obtained based on the received pilot symbols. 
[1050] The probabilities for the joint distribution Pr (0 = 0^,1 = Iq\p) may be 

expressed in the log-domain to simplify subsequent computation. This is similar to the 
use of log-likelihood ratio (LLR) to express the probability distribution of a single code 
bit. The use of log-domain representation for probabilities avoids the need to compute 
the outer exponential in the third expression of equation set (3). 

[1051] The joint probability distribution on 0 and / is obtained based on two 
variables 0 and / that are a priori independent. The joint distribution is thus a product 
distribution of a distribution on ^and a distribution on /. The distribution on 0 may be 
assumed to be uniform. The distribution on / (i.e., the interference power distribution) 
may be derived in various manners. In one embodiment, the interference power is 
assumed to be uniform over Q values. In another embodiment, the interference power 
is assumed to have a standard distribution such as a Gaussian distribution or a log- 
normal distribution. In yet another embodiment, the interference power distribution is 
obtained based on network-level computer simulation, empirical measurements, or by 
other means. 

[1052] In step 520, the forward LLRs for the B code bits of each received data 
symbol are computed based on the received data symbol and the joint probability 
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distribution on ^and / for that symbol. Step 520 is performed by each detection node 
420 in FIG. 4. The forward LLR computation for each detection node 420 may be 
performed in two steps. In the first step, an a posteriori distribution on the value of the 
transmitted data symbol, jc, is computed based on the received data symbol r and the 
joint probability distribution on 6 and /. A data symbol can take on one of M = 2® 
possible values. This a posteriori distribution indicates the probability for each of the 
M possible values of x given the received data symbol r and the joint probability 
distribution on dand /. The a posteriori distribution on x given r may be expressed as: 

Vr{x = x^\r) cx:Pr(r|^ = Xo) ^ Eq (4) 

where x is the data symbol value and can take on M possible values; 

Pr (jc = I r) is the a posteriori distribution on x given r, which gives the 

probability of the data symbol value being x^ given the received data 
symbol r; 

Pr (r I = jCq) is a distribution that gives the probability of obtaining received 

data symbol r given the data symbol value being x^ ; 

Pt{0 = 0^,I = Ij) is the joint distribution on 0 and / from the channel and 

interference estimation; and 
Pr (r I X = jCq,^ = d.,I = I j) is a complex Gaussian distribution with a mean of 

x^ae^^' and variance of Ij . 

The a posteriori distribution on x given r, Pr = j r) , may be viewed as a 2-D graph 

that contains M probability values for M values of x corresponding to the received data 
symbol r. 

[1053] In the second step for computing the forward LLRs, the a posteriori^ 
distribution on jc for the received data symbol is "marginalized" to obtain the forward 
LLRs for the B code bits of that symbol. The data symbol value x is determined by the 
values of the B code bits and the signal constellation used to map code bits to data 
symbols. Each of the M possible values of x corresponds to a different combination of 
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values for the B code bits. The forward LLR for each code bit can be computed as a 
weighted sum of the M probability values for the M possible values of x, where the 
weighting is determined by the distance between the data symbol in the signal 
constellation (or constellation symbol) and the received data symbol. Step 520 provides 
B forward LLRs for the B code bits of each received data symbol, or a total of Nb 
forward LLRs for the Nb code bits of the No data symbols being processed iteratively. 
The forward LLRs represent forward information provided to the decoder. 
[1054] In step 530, the forward LLRs for a codeword that contains the Nb code bits 
are decoded to obtain feedback LLRs for the Nb code bits. Step 530 is performed by 
decoder 340 in FIG. 4. The decoding may be performed based on, for example, a MAP, 
SOV, or Turbo decoding algorithm, and may be performed for one or multiple 
iterations. Step 530 provides Nb feedback LLRs for the Nb code bits, which represent 
feedback information provided by the decoder. 

[1055] In step 540, an a posteriori joint probability distribution on 0 and / is 
computed for each received data symbol based on the feedback LLRs for the B code 
bits of that data symbol, as follows: 

Pr (0 = 0^ J = I^\r) ocPr(/ = /,).Pr(r|^ = ^„/ = /,) , Eq (5) 

M 

= Pr (/ = /o) Y^Pr (X = X,)- Pr (r I X = x„0= 0„I = I^) , 

1=1 

where Pr (0 = 0q,I = I^lr) is the a posteriori joint probability distribution on 0and /, 

which provides the probability of the channel phase being 0^ and the 

interference power being 1^ given the received data symbol r; 

Pr(x = jc,.) is a product distribution on the data symbol value x, which can be 

obtained from the feedback LLRs for the B code bits of the received data 
symbol r; and 

Pr(r\x = x-,0 = 0Q,I = Iq) is a corhplex Gaussian random distribution with a 
mean of x- a- e^^° and variance of /q . 

Step 540 is performed by each estimation node 440 in FIG. 4. The a posteriori joint 
probability distribution on 0 and / given r is similar to the prior joint probability 
distribution on ^and / given p computed in step 510. However, the a posteriori joint 
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probability distribution is computed based on the feedback LLRs whereas the prior joint 
probability distribution is computed based on the received pilot symbol. Step 540 
provides No joint probability distributions on ^and / for No received data symbols. 
[1056] In step 550, the various joint probability distributions on 6 and / are 
combined to obtain updated joint probability distributions on ^and / for the No received 
data symbols. Step 550 is performed by node 410 in FIG. 4. For step 550, Nd a 
posteriori joint probability distributions on ^and / for the Nd received data symbols are 
available from step 540 and one prior joint probability distribution on 0 and / is 
available from step 510. These Nj, +1 distributions on 0 and / are used to derive Nd 
updated distributions on 0 and / for the Nd received data symbols. To avoid positive 
feedback, only extrinsic information is used to derive the updated distribution on 03x\d I 
for each detection node 420. The extrinsic information for a data symbol excludes 
information derived based on that data symbol. The updated distribution on ^and / for 
each received data symbol r is thus derived based on (1) Nj^-l a posteriori 
distributions on ^and / obtained for the other -1 received data symbols and (2) the 
prior distribution on 0 and / obtained from the received pilot symbols. This 
computation (1) effectively replaces the a posteriori distribution on 0and I obtained for 
received data symbol r with a uniform distribution and (2) assumes that the N^ - la 
posteriori distributions on 0 and / for the other Nj^— 1 received data symbols are 
obtained based on independent pieces of information. 

[1057] As an example, suppose that two distributions on 0 and /, namely 
Pr (0 = ^0,/ = /(, I a) and Pr (0 = 0^,1 = Iq\J5) , are to be combined, where aand are 

independent random variables when conditioned on 0and I. The composite distribution 
may be expressed as: 

Pr (0 = 0^ J = I^\a,j3) ocPr(a,j3\0 = 0^,I = I^) 

= Pria\0 = 0,J = I,)'Pr(j3\0 = 0,J = I,) Eq (6) 
ocPT(0 = 0^J = I^\a) FT(0 = 0^J = I^\m . 

The computation in equation (6) may be extended so that any number of distributions on 
^and / can be combined. Step 550 provides Nd updated joint probability distributions 
on ^and / for the Nd received data symbols, which are used by detection nodes 420 to 
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update the forward LLRs in the next iteration. Step 550 concludes one complete 
iteration of the joint channel and interference estimation and decoding. 
[1058] In step 560, a determination is made whether or not to terminate the iterative 
channel gain and interference estimation and decoding. This decision may be made 
based on one or more termination criteria. For example, a termination criterion may be 
as simple as a predetermined number of iterations. If the answer is *no* for step 560, 
then the process returns to step 520 to update the forward LLRs for the code bits. 
Otherwise, the process terminates. Step 560 may be performed after step 530, so that 
steps 540 and 550 can be omitted for the last iteration. 

[1059] For simplicity, the computation of probability distribution is exphcitly 
shown only for equation (1) and is omitted for all other equations. The computation of 
the various probability distributions described herein is known in the art. 
[1060] The number of values taken by channel phase and interference determine the 
computational complexity of the various steps in FIG. 5. Since each joint probability 
distribution on 0 and / contains L Q probability values, the computation complexity is 
proportional to both L and Q, which are the number of quantized values for channel 
phase and interference, respectively. To maintain reasonable complexity, interference 
may be quantized to a low resolution with a small number of values. A technique to 
reduce the number of channel phase values is described below. Complexity can also be 
reduced by using channel phase and interference in other manners, as described below. 
[1061] FIG. 6 shows a flow diagram of a process 600 for performing iterative 
channel estimation and decoding with an interference estimate. For process 600, 
interference is estimated once and used in the iterative channel estimation and decoding 
process. The interference estimate is not iteratively updated in order to reduce 
complexity. 

[1062] Initially, an interference estimate and prior information for channel phase are 
obtained based on the received pilot symbols (step 610). For step 610, a joint 
probability distribution on d and / is first obtained based on the received pilot symbols, 
as described above for step 510. The joint distribution on ^and / is then marginalized 
into a probability distribution on 6 and a probability distribution on /. A value of / is 
then selected based on the distribution on / and used as the interference estimate Ust- 
The interference estimate Ust may be the largest value in the distribution on /, the / value 
that results in minimum mean square error for the distribution on /, and so on. Step 610 
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provides a distribution on 6 and the interference estimate lesv The distribution on 6^ 
Pr = ^0 I P) ' "^^y expressed as: 

Pr(^ = ^olP)-Pr(p|^ = ^o./ = ^.J • Eq(7) 

Multiple distributions on 6 may be obtained for multiple received pilot symbols and 
combined to obtain one distribution on 0 for all received pilot symbols. Equation (7) 
represents one method of obtaining a distribution on d. Alternatively, the distribution 
resulting from marginalization can also be used directly to obtain the distribution on 
[1063] The forward LLRs for the B code bits of each received data symbol are then 
computed based on the received data symbol, the prior information for channel phase, 
and the interference estimate (step 620). Step 620 may be performed in two steps, 
similar to that described above for step 520 in FIG. 5. In the first step, an a posteriori 
distribution on x given r, Pr (jc = a:^ | r) , is computed as follows: 

Pr (x = I r) oc Pr (r I ;t = Xq) 

Eq (8) 

= £Pr(^ = ^.).Pr(r|x = Xo>^ = ^-^ = /eJ • 
1=1 

In the second step, the distribution on x for each received data symbol is marginalized to 
obtain the forward LLRs for the B codes bits of the received data symbol. 
[1064] The forward LLRs for all Nb code bits are then decoded to obtain feedback 
LLRs for the code bits (step 630). A posteriori information for channel phase is then 
obtained for each received data symbol based on the feedback LLRs for the B code bits 
of that symbol and the interference estimate (step 640). The a posteriori distribution on 
^for each received data symbol may be expressed as: 

Pr(^=^ok) -Pr(r|^ = ^o) 

Eq(9) 

»=i 

[1065] The a posteriori information for channel phase for the No received data 
symbols and the prior information for channel phase are then combined to obtain 
updated information for channel phase for each received data symbol (step 650). In 
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particular, an updated distribution on ^is computed for each received data symbol based 
on —la posteriori distributions on ^for the other —1 received data symbols 

and the prior distribution on ^derived from the received pilot symbols. 
[1066] A determination is then made whether or not to terminate the iterative 
channel estimation and decoding (step 660). If the answer is 'no', then the process 
returns to step 620 to update the forward LLRs for the B code bits of each received data 
symbol based on the received data symbol, the updated distribution on 0, and the 
interference estimate lest, as shown in equation (8). Otherwise, the process terminates. 
Again, step 660 may be performed after step 630. 

[1067] The iterative channel estimation and decoding with an interference estimate 
may also be performed in other manners, and this is within the scope of the invention. 
For example, the iterative process can start by computing a joint probability distribution 
on 0 and / for each received data symbol with the feedback LLRs for the symbol set to 
zeros. The received data symbol can provide information about ^modulo 2;r/M even 
if no information is available about the transmitted data symbols. This information 
about ^modulo Iniyi can then yield non-trivial information about the interference. 
The joint probability distributions on ^and / for the No received data symbols are then 
combined with a prior joint distribution on 0 and / derived from the received pilot 
symbols. The updated joint distribution on ^and / for each received data symbol is then 
marginalized to obtain a distribution on 6 and a distribution on /. A value of / is 
selected based on the distribution on / and used as the interference estimate lest^ The 
distribution on ^for each received data symbol and the interference estimate lest are then 
used to compute the forward LLRs for the received data symbol, as described above. 
[1068] As noted above, one factor that affects the computational complexity of a 
scheme that iteratively updates channel phase information is the number of values to 
which the channel phase ^is quantized (i.e., the value for L). For M-PSK modulation, 
the channel phase ^can be estimated to within a range of zero to 2niyi using non- 
iterative data-aided estimation. The technique described above for obtaining a 
distribution on ^for each received data symbol with the feedback LLRs set to zero is an 
example of data-aided estimation, which is a non-iterative technique because the code 
output was not used for estimation. The channel phase can then be quantized to M 
different values (instead of L values), which can greatly reduce computational 
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complexity if M is much less than L. Various methods may be used to estimate the 
channel phase based on the received data symbols. Two exemplary methods are 
described below. 

[1069] FIG. 7A illustrates the estimation of the channel phase Abased on the phases 
of the received data symbols. In this example, QPSK modulation is used, and five 
received data symbols are represented as solid dots 712a through 712e on a QPSK 
signal constellation 700. Each received data symbol has a phase that is determined by 
the modulation on the data symbol, the channel phase 0, and noise and interference. For 
this method, the phase of each received data symbol is first determined. A modulo- 
Iniyi operation is then performed on the phase of each received data symbol to 
effectively remove the modulation on the data symbol and convert the data symbol to a 
pilot symbol. The phases modulo- 2;?r/ M for the five received data symbols are plotted 
as circles 722a through 722e in signal constellation 700. The phases modulo- 2^/ M of 
the five data symbols are then averaged to obtain a channel phase estimate Oesty which is 
denoted by a "x" mark 730. 

[1070] FIG. 7B illustrates the estimation of the channel phase based on the complex 
values of the received data symbols. For this method, the five received data symbols 
are also represented as solid dots 712a through 712e on a QPSK signal constellation 
750. Each received data symbol is rotated by an integer multiple of Iniyi (i.e., by 
271 'i/M, where / is an integer zero or greater) so that the phase of the rotated received 
data symbol lies within a range of zero to 2^/M. Received data symbol 712a is 
rotated by zero degrees. Received data symbols 712b through 712e are rotated by the 
proper amounts and represented by circles 752b through 752e, respectively, in signal 
constellation 750. The five rotated received data symbols are averaged to obtain an 
average received data symbol, which is represented by a "x" mark 760. The channel 
phase 0 is estimated as the phase of the average received data symbol. This method 
uses both the amplitude and phase of the received data symbols to compute the channel 
phase estimate 0est, whereas the method illustrated in FIG. 7 A uses only the phase of the 
received data symbols. 

[1071] The channel phase may also be estimated in other manners, and this is within 
the scope of the invention. For example, the channel phase may be estimated using the 
techniques described above for obtaining a distribution on 0 for each received data 
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symbol with the feedback LLRs set to zero. As another example, phase of each 
received data symbol may be multiplied by M and then average. 

[1072] If the channel phase modulo- 2^ / M can be estimated reliably, then channel 
phase is known to within a range of zero to InIM, The only uncertainty is which one 
of M possible ranges the channel phase belongs to. An iterative scheme can then 
operate on M different values for channel phase, i.e., { 6^^^, 0^^^ +0^ , 0^^^ +2^^ , ... 

0^^^ -\- (M - 1) ' 0^ }, where - Itt/M, instead of L different values, where L is 
typically much larger than M. A distribution on ^ would then contain only M values or 
components, and a distribution on 0 and / would contain only M Q components. The 
computations for the forward LLRs and feedback LLRs are greatly simplified because 
of the fewer number of components to evaluate. For example, if the channel phase is 
quantized to L = 8M values, then the estimation of the channel phase to modulo- 

Iniyi reduces the channel estimation complexity by a factor of eight. The estimation 
of the channel phase (e.g., with received data symbols) may be used for any iterative 
receiver processing scheme. 

[1073] FIG. 8 shows a flow diagram of a process 800 for performing iterative 
channel estimation and decoding with a channel phase estimate obtained from received 
data symbols. For simplicity, process 800 does not account for interference. 
[1074] Initially, the channel phase ^is estimated to within a range of zero to 27r/ M 
based on the received data symbols (step 808). Step 808 may be performed with the 
first or second channel phase estimation method described above and provides the 
channel phase estimate . Prior information for channel phase is then obtained based 
on the received pilot symbols for M (instead of L) different values of ^(step 810). Step 
810 provides a distribution on Abased on the received pilot symbols for M different 
values of 0. 

[1075] The forward LLRs for the B code bits of each received data symbol are then 
computed based on the received data symbol and the prior information for channel 
phase (step 820). Step 820 may be performed in similar manner as described above for 
step 520 in FIG. 5 or step 620 in FIG. 6. The forward LLRs for all Nb code bits are 
decoded to obtain feedback LLRs for the code bits (step 830). A posteriori information 
for channel phase is then obtained for each received data symbol based on the feedback 
LLRs for the B code bits of that symbol and for M different values of ^(step 840). The 
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a posteriori information for channel phase for No received data symbols and the prior 
information for channel phase are combined to obtain updated information for channel 
phase for each received data symbol, again for M different values of ^(step 850). 
[1076] A determination is then made whether or not to terminate the iterative 
channel estimation and decoding (step 860). If the answer is 'no', then the process 
returns to step 820 to update the forward LLRs for the B code bits of each received data 
symbol based on the received data symbol and the updated distribution on 0. 
Otherwise, the process terminates. Again, step 860 may be performed after step 830. 
[1077] FIGS. 5, 6 and 8 show three specific schemes for performing iterative 
receiver processing. The scheme in FIG. 5 iteratively updates information for both 
channel phase and interference, the scheme in FIG. 6 iteratively updates information for 
channel phase and non-iteratively estimates interference, and the scheme in FIG. 8 
iteratively updates information for channel phase and uses phase information obtained 
from the received data symbols. Various other iterative schemes may also be 
implemented, and this is within the scope of the invention. For example, information 
for both channel phase and interference may be iteratively updated for a small number 
of iterations, and information for only channel phase may be iteratively updated 
thereafter. As another example, information for channel phase may be obtained once, 
and information for interference may be iteratively updated. As yet another example, 
the number of values for channel phase and interference can be reduced in subsequent 
iterations. As the distribution becomes more compact, some of the points have 
negligibly low probabilities and can be ignored. 

[1078] Referring back to FIGS. 3 and 4, the iterative channel and interference 
estimation and decoding scheme may be viewed as iterating between decoder 340 and 
an a posteriori probability (APP) detector 360 (with intervening channel deinterleaver 
336 and channel interleaver 342). (An APP detector returns a distribution on a bit (i.e., 
an LLR) whereas a MAP detector returns the most likely value of the bit (i.e., 0 or 1). 
APP detector 360 uses the feedback information (incoming LLRs) from decoder 340 
and the received pilot and data symbols (received values) to obtain the forward 
information (updated LLRs) for decoder 340. APP detector 360 estimates the channel 
and interference based on the feedback information and received values, and the channel 
and interference information is reflected in the forward information provided by APP 
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detector 360 to decoder 340. Other types of detector known in the art may also be used 
for the APP detector. 

[1079] In Tanner graph 400, channel and interference estimate node 410, detection 
nodes 420, and estimation nodes 440 represent one implementation of MAP detector 
360. MAP detector 360 may also be implemented in other manners, and this is within 
the scope of the invention. 

[1080] Referring back to FIG. 3, RX data processor 330 may implement any 
receiver processing scheme. For the scheme shown in FIG. 5, estimator 332 performs 
steps 510, 540 and 550, detector 334 performs step 520, and decoder 340 performs step 
530. For the scheme shown in FIG. 6, estimator 332 performs steps 610, 640 and 650, 
detector 334 performs step 620, and decoder 340 performs step 630. For the scheme 
shown in FIG. 8, estimator 332 performs steps 808, 810, 840 and 850, detector 334 
performs step 820, and decoder 340 performs step 830. 

[1081] For clarity, the processing by RX data processor 330 for the scheme shown 
in FIGS. 4 and 5 is described below. Channel and interference estimator 332 
implements node 410 and estimation nodes 440a through 440n in FIG. 4. For the first 
iteration, channel and interference estimator 332 obtains received pilot symbols from 
OFDM demodulator 320, derives prior information for channel phase and interference 
based on the received pilot symbols, and provides the prior information to detector 334. 
For each subsequent iteration, estimator 332 obtains received data symbols from OFDM 
demodulator 320 and feedback LLRs for the B code bits of each received data symbol 
from channel interleaver 342, derives a posteriori information for channel gain and 
interference for each received data symbol, combines the a posteriori information and 
the prior information, and provides updated information for channel gain and 
interference for each received data symbol to detector 334. 

[1082] Detector 334 implements detector nodes 420a through 420n in FIG. 4. 
Detector 334 obtains the received data symbols from OFDM demodulator 320 and the 
prior/updated information for channel phase and interference from estimator 332. 
Detector 334 computes the forward LLRs for the B code bits of each received data 
symbol based on the prior/updated information and the received data symbol and 
provides the forward LLRs to channel deinterleaver 336. Channel deinterleaver 336 
deinterleaves the forward LLRs. Decoder 340 performs decoding on the deinterleaved 
forward LLRs and provides feedback LLRs for the B code bits of each received data 
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symbol to channel interleaver 342. Channel interleaver 342 interleaves the feedback 
LLRs and provides interleaved feedback LLRs to estimator 332. 

[1083] In the description above, information for channel phase and interference is 
represented with probability distributions. Furthermore, forward and feedback 
information is represented with LLRs. Other representations may also be used for 
channel phase and interference and code bits, and this is within the scope of the 
invention. For example, log domain representation, inverse probabilities, and so on may 
be used. The computation for the various steps in FIGS. 5, 6, and 8 is dependent on the 
particular representations used for the channel phase and interference and code bits. 
[1084] For clarity, the computation for each of the steps in FIG. 5 is specifically 
described. Each step may also derive its information in other manners, and this is 
within the scope of the invention. As an example, for step 510 in FIG. 5, the received 
pilot symbols may be processed to obtain a pilot estimate, which is indicative of 
normalized channel observation for the pilot symbols, as shown in equation (2). The 
pilot estimate may be obtained by (1) multiplying the received pilot symbols with the 
conjugated pilot symbol to remove pilot modulation and (2) maximally combining all 
pilot symbols received for the symbol block being recovered. The pilot estimate is then 
used to obtain a joint probability distribution on 6 and / for the prior information for 
channel phase and interference. 

[1085] For simplicity, the description above assumes that the channel magnitude is 
known by the receiver (i.e., determined by some means) and only the channel phase is 
estimated iteratively. The channel magnitude may also be estimated iteratively along 
with the channel phase, albeit with an increase in complexity. 

[1086] The iterative receiver processing techniques described herein may be used 
for various wireless communication systems such as an OFDM-based system, a 
multiple-input multiple-output (MIMO) system, and so on. These techniques may also 
be used for the downlink (i.e., forward link) and the uplink (i.e., reverse link). For the 
downlink, transmitter 200 is part of an access point or a base station, and receiver 300 is 
part of a user terminal or a remote station. For the uplink, transmitter 200 is part of a 
user terminal, and receiver 300 is part of an access point. 

[1087] The iterative receiver processing techniques described herein may be 
implemented by various means. For example, these techniques may be implemented in 
hardware, software, or a combination thereof. For a hardware implementation, the 
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processing units used to perform iterative receiver processing (e.g., RX data processor 
330 in FIG. 3) may be implemented within one or more application specific integrated 
circuits (ASICs), digital signal processors (DSPs), digital signal processing devices 
(DSPDs), programmable logic devices (PLDs), field programmable gate arrays 
(FPGAs), processors, controllers, micro-controllers, microprocessors, other electronic 
units designed to perform the functions described herein, or a combination thereof. 
[1088] For a software implementation, the iterative receiver processing techniques 
may be implemented with modules (e.g., procedures, functions, and so on) that perform 
the functions described herein. The software codes may be stored in a memory unit 
(e.g., memory unit 352 in FIG. 3) and executed by a processor (e.g., controller 350). 
The memory unit may be implemented within the processor or external to the processor. 
[1089] The previous description of the disclosed embodiments is provided to enable 
any person skilled in the art to make or use the present invention. Various 
modifications to these embodiments will be readily apparent to those skilled in the art, 
and the generic principles defined herein may be applied to other embodiments without 
departing from the spirit or scope of the invention. Thus, the present invention is not 
intended to be limited to the embodiments shown herein but is to be accorded the widest 
scope consistent with the principles and novel features disclosed herein. 



[1090] WHAT IS CLAIMED IS: 



